Quantum criticality of spin-1 bosons in a ID harmonic trap 
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We investigate universal thermodynamics and quantum criticality of spin-1 bosons with strongly 
repulsive density-density and antiferromagnetic spin-exchange interactions in a one-dimensional har- 
monic trap. From the equation of state, we find that a partially-polarized core is surrounded by two 
wings composed of either spin-singlet pairs or a fully spin-aligned Tonks-Girardeau gas depending on 
the polarization. We describe how the scaling behaviour of density profiles can reveal the universal 
nature of quantum criticality and map out the quantum phase diagram. We further show that at 
quantum criticality the dynamical critical exponent z = 2 and correlation length exponent v — 1/2. 
This reveals a subtle resemblance to the physics of the spin- 1/2 attractive Fermi gas. 
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Alkali bosons with hyperfine spins in an optical trap 
provide exciting opportunities to simulate a variety of 
macroscopic quantum phenomena. In the spinor gas, the 
spin-dipolar collisions significantly change the spin states 
producing rich Zeeman effects. In particular, spinor 
Bose gases with density-density interaction and antiferro- 
magnetic spin-exchange interaction (2|] exhibit various 
phases of strongly correlated quantum liquids and are 
thus particularly valuable to study quantum magnetism 
and criticality. The experimental study of quantum crit- 
icality and universal scaling behaviour has recently been 
initiated in low-dimensional cold atomic matter 0, 
These advances build on theoretical schemes for map- 
ping out quantum criticality in cold atom systems [5|-l7|. 
In this framework, exactly solved models of cold atoms, 
exhibiting quantum phase transitions, provide a rigorous 
way to investigate quantum criticality Q . 

One-dimensional (ID) spinor Bose gases with short- 
range delta-function interaction and antiferromagnetic 
spin-exchange interaction are particularly interesting due 
to the existence of various phases of quantum liquids as- 
sociated with exact Bethe ansatz solutions 
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antiferromagnetic interaction leads to an effective attrac- 
tion in the spin-singlet channel that gives rise to a quasi- 
condensate of singlet bosonic pairs when the external 
field is less than a lower critical field at zero tempera- 
ture. In this phase, the low energy physics can be char- 
acterized by a spin-charge separation theory of the U{1) 
Tomonaga-Luttinger liquid (TLL) describing the charge 
sector and a 0(3) non- linear sigma model describing the 
spin sector [T^j. However, if the external field exceeds 
an upper critical field, we have a solely ferromagnetic 
quasi-condensate of unpaired bosons with spins aligned 
along the external field. For an intermediate magnetic 
field, the spin-singlet pairs and spin-aligned bosons form 
a two-component TLL with a magnetization, see Fig. [T] 

Quantum critical phenomena in the spin-1 Bose gas as- 
sociated with phase transitions at zero temperature can 




FIG. 1. (Color online) Phase diagram of the spin-1 Bose gas 
at T = in terms of the rescaled axial radius R/a^N^^'^ versus 
the polarization P = Ni/N for N/{alc^) = 0.05. The lines 
_Ri and R2 denote the phase boundaries of vanishing density 
of singlet-pairs and vanishing density of spin-aligned bosons. 
The model exhibits three quantum phases: spin-singlet paired 
bosons S, ferromagnetic spin-aligned bosons F, and a mixed 
phase of the pairs and unpaired bosons. V stands for the 
vacuum. 



be explored by varying the magnetic field and chemical 
potential. In this paper we use the exact thermodynamic 
Bethe ansatz (TBA) solution for spin-1 bosons to illus- 
trate the microscopic origin of the quantum criticality of 
different spin states in a ID harmonic trap. 

The model and equation of state - We consider N par- 
ticles of mass m confined in ID to a length L with 5- 
interacting type density-density and spin-exchange inter- 
actions between two atoms. The Hamiltonian is [HQ 
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where 5*^ is a spin-1 operator with z-component (s = 
1, 0,-1), Ez = —HS^ is the Zeeman energy, with H the 
external field and the total spin in the z-component. 
The ID interaction parameters cg — {gQ+2g2) /i and C2 — 
(ff2 — 5o)/3 with gs — ~2}i? /mas- Here m is the particle 
mass and as represents the ID s-wave scattering length 
in the total spin S* = 0,2 channels. Following jlj, 14 1 
the efi'ective ID coupling constant gs can be expressed 
in terms of the known 3D scattering lengths. 

Extending Yang's method [3] , Cao et al. solved the 
model IH) for the case c = cq = C2. For convenience, we 
define the linear density n = N/L and use dimensionless 
quantities ji = n/ £b, h = H/eh, t = T /ei, andp = p/\c\e},- 
Here energy and length are measured in units of the bind- 
ing energy Sb = fi^c^/Wm and c~^, respectively. We set 
H = 2m = 1 in the following. 

The model ([T]) has two conserved quantities, and the 
total particle number N. The number of particles in a 
particular spin state (s — 1, 0, —1) is no longer conserved. 
For instance the scattering between two particles of spin 
s — ±1 can produce two particles of spin s ~ and vice 
versa. Therefore, the model can have two types of pairing 
between two spin ±1 atoms or between two spin-0 atoms, 
see [11, Sim. 

The ground state properties of the system were stud- 
ied in [10l4l2l |. Analytic results can be extracted from 



the TEA equations for strong coupling 7^1 and low 
temperature t ^ 1 (see for details). The pressure 
of the system can be written as p — Pi + 1)2 in terms 
of the pressure pi for unpaired and p2 for paired bosons 
(Boltzmann constant /ss = 1) where 
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Here = Li„(— e'^^/*) with j — 1,2 denotes the standard 
poly log function Lin{x), with 
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The terms /, = te-''^^e-^P^/^Io{2pi/t) with Io{z) = 

X^fc^o ■^^^S — extracted from the spin wave bound 
states (16i] . 

The above result for the pressure serves as the equation 
of state, from which thermodynamic quantities such as 
the particle density n, the density of unpaired bosons ni, 
the density of spin-singlet pairs n2 , the entropy, and the 
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FIG. 2. (Color online) The density distribution profiles for 
trapped bosons at zero temperature with N/a^c? ~ 0.05 for 
(a) P « 12.22% (or h = 0.49) and (b) P « 22.33% (or h = 
0.51). 



compressibility of the system can be determined through 
the standard thermodynamic relations. Moreover, the 
equation of state allows an exact description of TLLs 
and quantum criticality of the system. 

Phase diagram and density profiles.- In order to study 
quantum criticality in the spin-1 Bose gas in a ID har- 
monic trap, the equation of state ^ can be reformu- 
lated within the local density approximation (LDA) by 
a replacement ^{x) — /i(0) — ^tow^o;^ in which x is the 
position and is the trapping frequency, see [l3| • Using 
the dimensionless chemical potential jl{x) = /i(0) — 8i^, 
where x = x/a^\c\ is a rescaled distance and = 
yjh/mujx is the harmonic characteristic radius, the to- 
tal particle number and the polarization are given by 
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We extract the phase boundaries from the equation of 
state (ID) within the LDA (g]) in the hmit t -> 0, see Fig. 
[T] The line Ri {R2) indicates the vanishing density of un- 
paired bosons (spin-singlet pairs). The intersection of the 
phase boundaries occurs a,t h = 0.5 and ^ = 0.5 yield- 
ing the critical polarization Pc ~ 16.97%. This phase 
diagram resembles that of the spin- 1/2 attractive Fermi 
gas, recently confirmed by Liao et al. [l3|. Here the 
critical properties of the spin-1 Bose gas are described 
by two Tonks-Girardeau gases with masses m and 2m, 
which are likely to mimic the free fermions and bosonic 
bound pairs in the spin-1/2 attractive Fermi gas. 

Fig. [2] shows the density distribution profiles in the 
limit t — > for N j c?^(? ~ 0.05 for (a) low polariza- 
tion {P < Pc) P ~ 12.22% and (b) large polarization 
(P > Pc) P ~ 22.33%. We find a partially polarized 
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FIG. 3. (Color online) Quantum criticality for low polariza- 
tion. The rescaled density profiles are plotted within LDA 
near the critical points (a) fid and (b) /ic2 at different tem- 
peratures t = 0.005, O.Of, 0.015 and 0.02. All curves in- 
tersect at the critical points. The rescaled density vs (c) 
(/i(a;) - fJ.c4)/T'"' and (d) {fi{x) - ^,c2)/T''"' at different tem- 
peratures t = 0.005, O.Of, 0.0f5 and 0.02. All data points 
collapse onto a single curve. 



core surrounded by two wings composed of eitfier spin- 
singlet pairs or a fully spin-aligned Tonks- Girardeau gas. 
The threshold values of the phase boundaries are con- 
sistent with the phase diagram FiglTJ Our result con- 
vincingly indicates that spinor Bose gases in a ID har- 
monic trap constitute an excellent quantum simulator of 
strongly correlated TLLs which could be probed b y i n 
situ imaging, analogously to the trapped Fermi gas [18 1. 

Quantum criticality - The equation of state ^ , captur- 
ing the entire thermal fluctuation in the quantum critical 
regimes, illustrates the microscopic origin of the quantum 
criticality. The thermodynamic functions can be cast in 
a universal scaling form ^_9„ 2^ For example, the density 
and compressibility are given by 

n(/x,T,x)=no-fr^+i-^g(^^M^^, (5) 

with dimensionality d — 1. Here we find the scaling func- 
tions Q{x) — AaLii/2(— e^) and J^[x) = A^Li_i/2(— e^), 
from which follow the dynamical critical exponent z = 2 
and correlation length exponent v = 1/2 for different 
phases of the spin states. In the above, no, kq, Aq, and 
A/3 are constants associated with the background near 
the critical points (see details in [16|). Moreover, quan- 
tum criticality driven by the external field H leads to 
similar scaling, where now Q{x) = AgLii/2(— e^"^) and 
F{x) = A/Li_i/2(— e'^"'^) with a dimensionless factor Aq 
and the same critical exponents fl6'|. Here Ag./ are dif- 
ferent constants associated with critical points. 



The critical exponents appearing in the scaling func- 
tions are analytically determined from the exact Bethe 
ansatz solution. We observe that the ID spin-1 Bose 
gas belongs to the same universality class as spin-1/2 at- 
tractive fermions [1] due to the hard-core nature of the 
two coupled Tonks-Girardeau gases. This demonstrates 
a significant feature of ID many-body physics - quan- 
tum criticality involves a universal crossover from a TLL 
with linear dispersion to free fermions with a quadratic 
dispersion in the low energy physics 2l|. In addition. 



the critical exponents of ID spin-1 bosons have the same 
values as the critical exponents of the Mott transition in 
ID and 2D Bose-Hubbard models. 



see 
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The universal scaling behaviour of the homogenous 
system can be mapped out through the density profiles 
of the trapped gas at finite temperatures. Although a 
finite-size scaling effect is evident in the scaling analysis 
[S S US] , the finite-size error lies within the current ex- 
perimental accuracy [l^ . Here we show that all thermo- 
dynamic observables are given in terms of the rescaled 
temperature and chemical potential. Thus one can ei- 
ther lower the temperature or increase the interaction 
strength (by tuning the transverse frequency) such that 
all data curves for the physical properties in the trapped 
gas collapse at different temperatures into a single curve 
with a proper scaling. This collapse signature can be used 
to confirm the quantum critical law, see the experimental 
measurements given in [23| . 

In the regime of low polarization [P < Pc), the phase 
transitions from the vacuum into the spin-singlet paired 
phase and from the pure paired phase into the mix- 
ture of spin-singlet pairs and spin-aligned bosons oc- 
cur as the chemical potential passes the lower critical 
point jlc2 = —0.5 and the upper critical point fid — 
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The universal scaling behaviour can be identified 
following the scheme proposed in Q by plotting the 
"rescaled density" versus the chemical potential for dif- 
ferent values of temperature. Figures [3] (a) and (b) show 
the rescaled density normalized by N-^/^/{a,j;) versus the 
rescaled distance at different temperatures. All curves 
intersect at the critical points /ic2 and ftci after a sub- 
traction of the background density no. We also see in 
Fig. [3] that all data points collapse onto a single curve 
near (c) /ic4 and (d) /ic2 , which is the scaling function for 
these quantum critical points. 

For large polarization (P > P^), the phase transi- 
tions from vacuum into the ferromagnetic spin-aligned 
boson phase and from the spin-aligned boson phase into 
the mixture of spin-singlet pairs and spin-aligned bosons 
occur as the chemical potential varies across the lower 
critical point p,ci = —h and the upper critical point 

Mc3 = +'^ {h-^S' +7B {h-iy> respectively. 
A similar analysis is presented in Fig. 2] near the critical 
points fici and jlcs- 
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FIG. 4. (Color online) Quantum criticality for high polariza- 
tion. The rescaled density profiles are plotted within LDA 
near the critical points (a) fid and (b) /ic3- The rescaled den- 
sity vs (c) {fi{x) - Mci)/y and (d) - /ic3)/y • As 
in Fig. the curves and data points are indicative of scaling 
behaviour at quantum criticality. 



Similar plots can be made for the rescaled compress- 
ibility in agreement with the universal scaling result ([6]) . 
The compressibility curves at different temperatures col- 
lapse into a single curve obeying the scaling law. The 
compressibility always develops a peak near the phase 
boundary on the side with higher density of state. How- 
ever, the experimental measurement of the compressibil- 
ity is more difficult than the density profiles. Neverthe- 
less, progress in measuring the compressibility of an ul- 
tracold Fermi gas has recently been made 25 1. 

In conclusion, we have discussed the equation of state, 
the density profiles, TLL thermodynamics and univer- 
sal scaling at quantum criticality for the spin-1 Bose gas 
with strongly repulsive density-density and antiferromag- 
netic spin-exchange interactions in a f D harmonic trap. 
We have shown that the quantum criticality of different 
spin states belongs to the universality class with criti- 
cal exponents z = 2 and v = 1/2. This remains true 
in general for the ID model even for the non-integrable 
case Co 7^ C2. We have also demonstrated that the phase 
diagram, the TLLs and critical properties of the bulk 
system can be mapped out from the density profiles of 
the trapped spinorgas at finite temperatures. Current 
experiments [3, 0, E, [l^ are capable of measuring such 
universal features of ID many-body physics. 
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